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Abstract 


Differential cross sections for electromagnetic dissociation in nuclear collisions 
are calculated for the first time. In order to be useful for three - dimensional 
transport codes, these cross sections have been calculated in both the projectile 
and lab frames. The formulas for these cross sections are such that they can 
be immediately used in space radiation transport codes. Only a limited amount 
of data exists, but the comparison between theory and experiment is good. 


1 Introduction 

The issue of protecting astronauts from cosmic radiation is becoming increasingly 
important with the current plans to establish a permanent human base on the Moon 
with a follow-on mission to Mars. During a Mars mission, astronauts will be exposed 
to cosmic radiation fields for several years. What is unique to cosmic radiation fields is 
the presence of heavy nuclei (up to and beyond Z=26) with energies of GeV and beyond. 
In calculating space radiation environments, it is therefore necessary to have a detailed 
knowledge of nucleus - nucleus collisions in the GeV energy region where the cosmic ray 
spectrum reaches a maximum. 

Nucleus - nucleus collisions are mediated mainly by either the strong or electromag- 
netic (EM) force. A reaction proceeding via the EM force is often called Electromagnetic 
Dissociation (EMD). For few nucleon removal, EMD cross sections are just as important 
as strong interaction cross sections. Differential EMD cross sections are necessary because 
fully three-dimensional transport codes require energy and angular differential cross sec- 
tions. A full theory of differential EMD cross sections has never been worked out. That 
is the goal of the present paper, which is an expanded version of Reference [1]. 


2 Photonuclear cross sections 

Nucleus - nucleus electromagnetic cross sections require photonuclear cross sections 
as input. Therefore, in this section, we first discuss these simpler and more basic pho- 
tonuclear processes. Included is a discussion of the differential photonuclear cross section, 
which will be important when determining the nucleus - nucleus differential EMD cross 
section. We also present a review of the total cross section parameterization, because this 
will be used in obtaining differential cross sections. 

2.1 Compound nucleus 

Some excellent references that discuss compound nucleus formation are Frobrich [2] 
(pp. 372, 322, 346), Preston [3] (pp. 514-515 ), Blatt and Weisskopf [4] (p. 368), Paul [5] 
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(pp. 220-224), Burcham and Jobes [6] (p. 189), de Benedetti [7], Bethe and Morrison [8], 
Harakeh and van der Woude [9], Satchler [10] (pp. 75, 249), Feshbach [11] (pp. 296, 309, 
326), McCarthy [12] and Jackson [13] (p. 94). 

An essential point emphasized by Burcham and Jobes [6] (p. 191) is that the compound 
nucleus is a true intermediate resonance state, so that the reaction proceeds according to 

X + a -> C* -> Y + b, (1) 

where C* represents the intermediate compound nucleus, rather than a direct reaction 

X + a -»• Y + b. (2) 

This means that the kinematics are entirely different. For the direct reaction (2), one 
cannot form a double differential cross section for the final particle b. Instead, || 
and ^ are not independent, but are functions of each other. However, for the compound 
nucleus reaction, an essential feature is that the compound nucleus “forgets” how it was 
formed after reaching statistical equilibrium. This means that it is a true intermediate 
state and one instead focuses on compound nucleus decay 

C* Y + b. (3) 

Here, the spectral and angular distributions will be independent of each other, so that 
one can now form the double differential cross section ^Edn ^ or ^ ie particle b. In the 
simplest model, the angular distributions will be constant and the spectral distribution 
will represent the statistical thermal decay of a Boltzmann system. The book by Preston 
[3] contains one of the most useful discussions of the compound nucleus, and a discussion 
of for the compound nucleus decay can be found on pages 514 - 515. 

2.2 Two - body final state 

We now consider reactions involving two or three particles in the final state. Suppose 
we have a reaction with only two bodies in the final state, such as 

1 + 2 -> 3 + 4 , (4) 

where the number denotes the particle. An example of such a reaction is 

7 + A — > A* — > {A — l)o + N , (5) 

where 7 denotes a photon, A denotes a parent nucleus, N denotes a nucleon, and ( A — 1) 0 
denotes a daughter nucleus in its ground state. In this reaction, the incident photon 
excites the parent nucleus to a compound nucleus excited state denoted by A*. The 
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excited parent compound nucleus decays directly to the ground state of the daughter, 
with the emission of a nucleon, as shown in figure 1. 


2.3 Three - body final state 

Now consider a three - body final state, such as 

1 + 2— >3 + 4 + 5, (6) 

where the number denotes the particle. An example of such a reaction is 

7 + A — > A* — > (A — 1)* + ./V — > (A — 1) 0 + 7 + ./V. (7) 

In this reaction, the parent nucleus is excited to an intermediate compound nucleus state, 
but now it decays to an excited state of the daughter nucleus again with the emission of a 
nucleon. The excited daughter then decays to its ground state by emitting a photon. (It 
could also decay to another excited state and emit another photon.) As shown in figure 
2, there will be a variety of excited daughter states that are possible to populate. 

For photonuclear reactions, one can measure both and for an emitted nucleon 

N. Equivalently, one can form dE ~ dUls , ■ If the reaction occurred via equation (5), then 
both differential cross sections cannot be determined. It is important to realize that 
both angular and spectral distributions are only possible when the photonuclear reaction 
involves a three - body final state as in equation (7). Alternatively, this reaction could be 
thought of as involving a two - body final state 7 + A — > (A — 1)* + N where the mass 
of the daughter nucleus (A — 1)* is not fixed but can vary, because it can exist in many 
different excited states or energy levels as shown in figure 2. This means effectively that 
the mass of the excited daughter nucleus is not fixed. Thus, the energy of the emitted 
nucleon can vary allowing the formation of 4+ 


2.4 Angular distribution 

In the simplest compound nucleus model, the photonuclear angular distribution is 
approximately isotropic [3, 5, 10, 14], meaning that the angular distribution is constant, 


djCJryA ( TV ) 


dQ 


= K 


N 


( 8 ) 


where K is a constant. It is trivial to evaluate K from the total cross section because 


ct 1 a(E 1 ) = f <lcr ~>^ E ^ d,Q N = AttK. 

J CLILn 


( 9 ) 
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Rearranging gives 


_ Cr 'yA(Ery) 

47T 

An isotropic angular distribution can be written as 

d(7~, A (Eyj (Ay A (Eyj 

dQ n 47 r 


( 10 ) 


( 11 ) 


2.5 Spectral distribution 

The energy level density in figure 2 can be approximated [2], [11] (p. 326) by a 
Boltzmann distribution 


p(E) ~ e ~ E / k& , (12) 

with the nuclear temperature given by [2, 11] 

“ = /?• < 13 > 

where d is a constant, sometimes taken as d— 10. Here, 0 is the nuclear temperature and 
k is the Boltzmann constant. The photonuclear spectral distribution is parameterized as 


d(JryA (.Ey) 
dE N 


C T n e~ TN/ke , 


(14) 


where Tn is the kinetic energy of the emitted nucleon. We write instead of 
because dE N = cITn ■ The total energy is given by En = Tn + tun, where tun is the 
nucleon mass. Also, note that we do not have i?jve~ Ejv / fc0 , because En begins at mjv 
rather than zero. 


2.5.1 Useful integrals 

In order to evaluate the constant C in the above expression for the spectral distribu- 
tion, we shall need the following integrals [15], [16] 


dx xe~ 


1 + ax 


J dE Ee~ E/ke = -k0(E + k0) e~ E/ie , 


(15) 

(16) 


4 



and 


/ dx x n e~ ax 
Jo 

r oo 

/ dx xe~ ax 

r dE Ee~ E/ke 
Jo 


T(n + 1) n\ 
a n+ 1 _ a n+1 ’ 

1 

9 ’ 

er 

(kef . 


( 17 ) 

(18) 
(19) 


However, we will need these integrals with different limits, such as a minimum or maximum 
energy, E min = -^threshold or E max . Thus, the following integrals will be useful. From 
reference [15], we obtain 


rC 


dx xe~ 


IB 


1 + aB 


1 + aB 


e~ aB - 


—aB 


1 + aC 


~aC 


for C = oo and Re [a] > 0. 


( 20 ) 


This reduces to the above result for B = 0 and C = oo, 


rB 'max ri/ir\ 

/ dE Ee~ E/ke 
JE min 


k<d{E min + ke)e~ E ^ /ke - kQ(E max + kQ)e~ E — /ke 
kO{E min + kG)e~ E ^ /ke , for £ max = oo. (21) 


2.5.2 Spectral distribution in terms of total cross section 

The photonu clear spectral distribution is written in equation (14), where Tn is the 
kinetic energy of the emitted nucleon, and C is a constant determined by the requirement 

-Ms,) = / d.E N . (22) 

where crt 0 t(-E 7 ) is the photonuclear total cross section. First, assume that the limits of 
integration are 0 and oo. Then, 

u tot (£ 7 ) = r dE N =C [°° dT N T N e~ T ^ ke = C(kG ) 2 , (23) 

Jo dE]y Jo 


giving 


O'tot ( E ~: ) 

(kQ) 2 


(24) 
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or 


da^A^Ery) _ q-tot(-^ 7 ) rp r N /ke 

dE N ~ (A:©) 2 N 


(25) 


which is the photonuclear spectral distribution. 

We see that this has the correct units, noting that the units of kQ are MeV. (In units 
where k = 1, 0 is in MeV. This is used in the computer codes for this work.) Equation 
(25) makes sense because the differential cross section is proportional to the total cross 
section times the probability of decay, ~ ae^ TN ^ ke . See figures 2 and 3 of reference 
[17]. The more general calculation, with arbitrary limits T min and T max , yields 


Otot (-^7) 



C 


k(~)(T„ 


d (T max 

— = C / dT N T N e- T »' kQ 
+ ke)e- T -~ /m - kO(T mm + /te)e- T "“*/‘ e l . 


(26) 


Rearranging gives 
C = 


Otot ( E~ ; 


kQ(T mm + kO)e~ T ^/ ke - A;0(T max + k<d)e~ T ™/ ke 
Substituting this value of C into equation (14) yields 
da 


(27) 


CTtot (.Ery 
- k( 

T n e~ TN/ke , 


dE N k(~) (T min + A;0)e-TninAe - kQ(T max + jfc0) e - W*e " N 

0" t ot {E- 


T n e ~ TN/ke 


for T max = 00 . 


A;0(T min + jfc0) e -W*e 
This reduces to the above result, equation (25), when T min = 0 and T max = 00 . 


(28) 


2.6 Double differential cross section 


From reference [14] (pp. 27, 40) (with / = 0), the photonuclear double differential 
cross section can be expressed as 


d 2 a 1 A(E, y ) 1 da^iEj) 

dEjsidElN 47t dE n 

which corresponds to an isotropic angular distribution. 


(29) 
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2.7 Lorentz invariant differential cross section 


The Lorentz invariant differential cross section E^ is related to the non-invariant 
double differential cross section via 


E (P<Ty A (Ej _ d?<T yA (Ej 
N d 3 p N p N (!E N dQ N 


(30) 


where p n = |pjvl- Note that the entire right hand side is to be evaluated in the same 
frame. For example, if the double differential cross section dE ^ dn is evaluated in the 
projectile frame, then the term also refers to the projectile frame. If the angular 
distribution is isotropic, we have 


(i (Jry A (Ey') 1 d(T ’ry A (^E iy ) 

d 3 pN 4 tt pn dE N 


(31) 


2.8 Total cross section 

The above equations for the photonuclear differential cross sections were all written 
in terms of the photonuclear total cross section. We now present relevant equations to 
calculate the photonuclear total cross section. The photonuclear total cross section for 
producing particle X is [18, 19] 


cr(E y ,X) = 3PfCr ab s(-E 7 ) , (32) 

where gx is the branching ratio and cr a b s (-E 7 ) is the photonuclear absorption cross section, 
which is parameterized as 


(7 abs ( E ~, ) 


(Jr 


1 + 


(£2 - e% dr )*/e* P 


(33) 


The abbreviation, GDR, stands for giant dipole resonance. Here, T/gdr. is the energy at 
which the photonuclear cross section has its peak value and T is the width of the electric 
dipole (El) giant dipole resonance. Also, 


(7 TRK 

with the Thomas-Reiche-Kuhn cross section given by [19] 


o'trk — 


QONpZp 

A P 


MeV mb , 


(34) 


(35) 
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with the subscript P referring to excitation of the projectile. (Note that in reference [19] 
a typing error had this subscript referring to the target.) The GDR energy is 


with 


Ecrm — 


he 


"m*c 2 i?Q 

8 J 


(1 + U 


1+e+u 



U = 


3 J 

Q 7 
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and 


fto = r 0 A^\ 


(36) 


(37) 


(38) 


(In reference [19] a typing error had this subscript referring to the target.) The parameters 
are: 


e = 

0.0768 , 

(39) 

Q' = 

17 MeV , 

(40) 

J = 

36.8 MeV , 

(41) 

r 0 = 

1.18 fm , 

(42) 

m* = 

0.7 m nuc leon • 

(43) 


3 Nucleus - nucleus cross sections 

In a nucleus - nucleus collision mediated by the EM force, the target (or projectile) 
represents a source of virtual photons, which impinge upon the projectile (or target). The 
spectrum of virtual photons contains a variety of energies, in contrast to a photonuclear 
reaction where the incoming photon possesses only a single energy. For a nucleus - nucleus 
collision, figures 1 and 2 get replaced by figures 3 and 4, which show a variety of energy 
levels being excited in the parent nucleus. 

3.1 Total cross section 

The total cross section for nucleus - nucleus reactions can be written in the form 

vaa = J dE 7 N(E 1 ) (t 7 a(E 7 ) , (44) 

where N(E^) is the Weiszacker- Williams virtual photon spectrum and <t 7 a(-E 7 ) is the 
photonuclear total cross section. 



3.2 Angular distribution 

From the perspective of EMD reactions, the spectator nucleus is nothing more than 
a source of virtual photons, therefore the angular and spectral distributions may also be 
written in the form of equation (44). Thus, the angular distribution, for emission of a 
nucleon N in the direction Qft, is given by 


dcr a A 


= J dE 7 N(E^) 


d.’CT 'yA ( -hf 


dQ 


N 


(45) 


where is the photonuclear angular distribution for emission of a nucleon N in the 

direction Qft. If the photonuclear angular distribution is approximately isotropic, then 
use of equations (11) and (45) gives 


d&AA _ &AA 
dElft 47T 


(46) 


It must be emphasized that this is the angular distribution in the rest frame of the com- 
pound nucleus. If the compound nucleus is the projectile, then this must be transformed 
to the lab frame for use in transport codes. If the compound nucleus is the target, then 
no transformation is necessary, because the target is at rest in the lab frame. 


3.3 Spectral distribution 

The spectral distribution, for emission of a nucleon N with energy Eft, may also be 
written in the form of equation (44), namely 

(4T) 

where — is the photonuclear spectral distribution for emission of a nucleon N with 
energy E N . Note that the spectral distribution cannot be taken outside the integral because 
the nuclear temperature 0 depends on the photon excitation energy E 1 . 


3.4 Double differential cross section 


The double differential cross section, for emission of a nucleon N with energy Eft in 
the direction Qft, may also be written in the form of equation (44), namely 


d 2 aAA 
dQ ftdEft 


= j dE 7 N(EJ 


d 2 a^ A (E y ) 

dQ ftdEft 


(48) 
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where is the photonuclear double differential cross section for emission of a 

nucleon N with energy E n in the direction If the photonuclear angular distribution 
is isotropic, we use equation (29) to give 


d 2 (T A A 

di'liKjdEpj 




1 d(JAA 

47 t d,E N 


(49) 


which is analogous to equation (29). 


3.5 Lorentz invariant differential cross section 


The Lorentz invariant differential cross section is related to the non-invariant double 
differential cross section as before, namely 

d 3 a AA 1 d 2 a a a 

— — jj? jo ' ( 50 ) 

d Pn Pn dE N d\l N 

For an isotropic distribution, we have the analog of equation (31), 


d 3 (T AA _ 1 da AA 

d 3 pN 4:7rpN dE N 


(51) 


4 Lorentz transformation of cross sections 

The nucleus - nucleus differential cross sections in equations (45), (47), (48) all involve 
a photonuclear differential cross section or a total cross section. The photonuclear dif- 
ferential cross sections are all evaluated in the rest frame of the nucleus undergoing the 
photonuclear reaction. Differential cross sections in radiation transport codes are required 
in the lab frame (spacecraft rest frame). If the projectile nucleus is undergoing photodis- 
integration, then the nucleus - nucleus differential cross sections in equations (45), (47), 
(48) are first evaluated in the projectile frame. Then, they must be Lorentz transformed 
to the lab frame. The technique for doing this is the subject of the present section. 


4.1 Discussion of photonuclear cross sections 

For a two - body final state, such as the reaction in equation (5), one can form either 
the angular or spectral distribution, but not both. A double differential cross section or 
Lorentz invariant differential cross section cannot be formed in this case. For a three - body 
final state, such as the reaction in equation (7), both the spectral and angular distributions 
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can be formed independently, as well as double and Lorentz invariant differential cross 
sections. This can be seen in the literature for photonuclear reactions. 

For general photonuclear reactions, numerous discussions of double differential cross 
sections [17], [20] - [27] can be found. However, for photonuclear reactions to the ground 
state, one can find information regarding only angular distributions [28] - [31]. This result 
is expected because when making a transition to the ground state, as shown in figure 1, 
there will be no energy distribution for the emitted nucleon when assuming the incident 
photon is monoenergetic. This is true for any transition to a particular energy state in the 
daughter nucleus. Although, when one allows for transitions to a variety of energy states 
in the daughter, as shown in figure 2, there will then be a variety of possible emitted 
nucleon energies whose spectral distribution will reflect the energy level density of the 
daughter nucleus. A transition to the ground state of the daughter nucleus is distinctive 
because no photon is emitted, and the final state will be a two - body state. Consequently, 
the angular distribution calculated in the projectile frame will need to be transformed into 
the lab frame, which is a complicated task. 

From the view of a statistical compound nucleus decay, there is nothing special about 
the ground state. It is just one of a continuum of possible final states, which is given 
by the continuous energy level distribution p(E). In this general case, the final state is 
three - body. This enables the formation of a Lorentz invariant differential cross section, 
which is easily transformed from the projectile to the lab frame. (Even though E ^ is 
invariant, it must be transformed when plotting it as a function of energy or angle.) 

4.2 Lorentz transformation of photonuclear cross sections 

A Lorentz transformation of a photonuclear differential cross section is almost never 
considered because it is inherently defined in the rest frame of the nucleus, which un- 
dergoes the reaction. The projectile is a photon, the target is a nucleus, and Lorentz 
transformation is not necessary. When considering nucleus - nucleus collisions, Lorentz 
transformations are needed. Assume that the projectile nucleus is undergoing the pho- 
tonuclear reaction. Now, we must transform the projectile photonuclear differential cross 
section into the target nucleus frame. 

In this section, when a Lorentz transformation of a differential photonuclear cross 
section is discussed, we are considering the different frames available to the nucleus - 
nucleus reaction, not the photonuclear reaction. This is admittedly a confusing point. It 
is important to realize that the photonuclear reaction is still occurring to the projectile 
nucleus. Normally, we would just integrate the projectile photonuclear cross section over 
the virtual photon spectrum coming from the target nucleus. We still do this, but first we 
transform the projectile photonuclear cross section to the lab (target) frame. It should be 
noted that although the photonuclear cross sections are transformed to the lab frame and 
written in terms of lab variables, they are still photonuclear cross sections for projectile 
fragmentation. 
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4.3 Lorentz transformation between cm or projectile frame and 
lab (target) frame 

Suppose we have quantities in the center of momentum (cm) or projectile frames and 
we wish to transform to the lab frame. The cm frame moves at speed (3 c i relative to the lab 
frame. The corresponding 7 factor is labeled as 7 d- The projectile frame moves at speed 
(3 p i relative to the lab frame. The corresponding 7 factor is labeled as 7 p /. The Lorentz 
transformations are 


( E* \ = ( 1 *1 -7 \ f Ei \ 

V p\\* ) V 1*1 ) V p\\i J ’ 

and inverse transformations are 


Pt* — Pti , 


where 



( l*i l*i(d*i \ ( E * \ 

\ l *i(3*i 1*1 ) \ P\\* ) ’ 


Pti — Pt* , 


P\\=Pz = |p|cos0, 
Pt = |p|sind. 


(52) 


(53) 


(54) 

(55) 


With this notation, both cm and projectile frames are included. The notation means that 
a quantity x * is the value of the quantity x evaluated in that particular frame with 


* = c or * = p, 


(56) 


where the c or p subscript refers to the cm or projectile frame, respectively, and /3*i is the 
speed of that frame with respect to the lab frame with 


(3*i = (3d or (3d = (3 p i . 


(57) 


4.4 Energy transformations 

The energy Lorentz transformation from the lab (/) frame to the starred (*) frame is 


Ej* id{Eji l*iP\ \ji ) 

= l*i(Eji - (3*i\pji\ cos Op) 

= 7 *) [Eji - (3*i\J E? t — m 2 cos 6 ^ 7 ) . (58) 

The inverse transformation is 

Eji = 7 */ (E jit + (3d\J T 7 * — rn 2 cos d,*) . (59) 
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4.5 Angle transformations 

The angle is obtained from 


tan 9 


Pr 

Pz 


(60) 


Thus, the angle of particle j is 

tan 9 ji = 

Pzjl 


Pt j* 

y*;/?*^* T 7*;p 2 j* 

|pj*| sin6(j* 

7 *i (/ %iEj * + |pj*| cos Oj*) 


(61) 


Defining <37* as the speed of the cm or projectile frame relative to the lab frame divided 
by the speed of particle j in the cm or projectile frame 

aj, = ?P, (62) 


and using 


flj* 


Pj*[ 

rp ’ 
^3* 


we obtain 


tan 6ji = 


sin Qj* 

7 *1 (cos 9 j* + 07*) 


(63) 


(64) 


See references [32] (p. 402), [33] (p. 42), [34] (p. 17), [35] (p. 26). This is a complicated 
function of 9 because, in general, 


a j* = a j*(EjP = 07* (0j*)- (65) 

Usually, 07* is a function of 9j*, making tan 9 3 i a complicated function of 9j*. For the cm 
frame and for three - body states, however, Ej C is not a function of 9j C . This means that 
atjc is not a function of 9j C [33] (pp. f2, 58). Three - body states are considered in the 
present work, so that we avoid the aforementioned complications. 

The angle transformation for * = c is plotted in references [32] (p. 403), [33] (p. 43), 
[36]). When 07 * > 1, the function is double valued. This results from the two different 
angles in the cm or projectile frame giving rise to the same angle in the lab frame for 
07* > 1. The two angles can be distinguished by their energies, labeled in the cm frame 
as Ef c [32] (p. 402). This is true for a two - body final state. 
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To show that the cm or projectile angle is double valued, we specify the * frame in 
equation (64) to be the projectile frame. Then, 


tan 9ji = 


sin 0j p 

'T pi (cos 9 jp T Qjp) 


( 66 ) 


where aj p is defined as the ratio of the projectile velocity to the velocity of particle j 
in the projectile system. It is dependent on the lab kinetic energy T/v, as in aq p (T/v). 
When aj P > 1, particles are emitted forwards and backwards in the projectile system, 
but appear at the same lab angle. Figures 5-8 are evaluated at the lab kinetic energy 
of mjv, lOm^v, 16mjv, and IOOttttv, respectively. Figures 5 and 6 are examples of what 
happens when aj P > 1. It can be seen that the range of the lab angle, 0 < 6i a b < 6 max , is 
dependent on a rp . The lab angle is also confined to a forward cone, as the projectile angle 
ranges from 0 to it. This means that the projectile angle will be double valued for every 
lab angle. Figures 7 and 8, however, illustrate the relationship between the projectile 
angle and the lab angle when a p < 1. Notice that the lab angle range is now between 0 
and 7 r. The lab angle is now a single valued function of the projectile angle. 

In the present work, we consider three - body final states. The curves plotted in 
reference [32] (p. 403) correspond to a particular value of the particle energy, or in other 
words, to a particular value of ay*. For a three - body reaction, the energy of the emitted 
particle is not related to the angle. Thus, for a particular lab angle, we have a continuous 
range of lab energies corresponding to a family of curves. This is plotted in reference [32] 
(p. 403). The family of curves will be for both a > 1 and a < 1. 

The inverse transformation is 


where 


- Sln fy 

71 ( 008 ^,- 0 ^) ’ 


pfl 


a tl = 


Ipjl 

E i 


*1 


ft 


(67) 


( 68 ) 

(69) 


The ± notation is emphasizing that, in general, two different angles in the cm or projectile 
frames can correspond to a single angle in the lab frame. 
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4.6 Double differential cross sections 


The transformation of a double differential cross section is described in reference [37] . 
The cross sections are related by the Jacobian 


which is evaluated as 


cPa d 2 cr d 

dEjidQji dEj^dflj* d (Eji, f Iji) 

d (i?*,n*) _ |p;| _ sin 6** 

d (Ei, £li) |p* | sin 9i 


(70) 


(71) 


to give 

d 2 a |p/| d 2 <j sin 67* d 2 a 

= | — - = - . (72) 

dEjidilji |p* | dEjzdQj* sin 9ji dEj+dilj* 

The left hand side is a function of Eji and dji, so that the right hand side also should 
be function of Eji and 6ji. This is accomplished by replacing Ej * and (7,-* on the right 
hand side with equations (58) and (67). Thus, the entire right hand side is written as an 
explicit function of lab variables. We have written two versions of the right hand side. 
Either version can be used, and the one involving the sine functions will be a good test 
of the correctness of the angle transformations. 

Equation (72) is the equation for obtaining all cross sections in the lab frame. The 
method is to use equation (72) to obtain the photonuclear double differential cross section 
in the lab frame. Then integrate the lab frame double differential cross section in equation 
(72), to get the photonuclear spectral and angular differential cross sections in the lab 
frame. To obtain any nucleus - nucleus differential cross section in the lab frame, take 
the lab frame photonuclear cross sections and integrate over the virtual photon spectrum. 


5 Results 

The results of these calculations are presented in figures 9 - 22. The photonuclear 
differential cross sections are presented in figures 9 - 14, which are evaluated at a photon 
energy of 20 MeV, near the peak of a typical giant resonance. The nucleus - nucleus 
differential cross sections are presented in figures 15 - 20. The cross sections in figures 9 
- 20 are for the reaction 


28 Si + 208 Pb -> n + 27 Si + 208 Pb (73) 

at 14.6 AGeV. Figure 9 shows the application of equation (25) to calculate the pho- 
tonuclear spectral distribution in the lab frame. Equation (11) was used to calculate an 
isotropic photonuclear angular distribution in the lab frame, as seen in figure 10. The 
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photonuclear double differential cross section in the lab frame is given by equation (29) 
and is shown in figure 11. 

These photonuclear cross sections are transformed to the nucleus - nucleus lab frame 
and are shown in figures 12 - 14. The photonuclear double differential cross section 
in the lab frame, as shown in figure 12, displays the double peak feature discussed by 
Hagedorn [38] (pp. 47 - 49). This double peak occurs because there is a single peak in the 
spectral distribution in the projectile frame (figure 9), which gets boosted both forward 
and backward in the lab frame depending on kinematic conditions. Figure 13 displays 
the photonuclear spectral distribution in the lab frame, which results from integrating 
the photonuclear double differential cross section in the lab frame over all lab angles. 
Compared to figure 9, it can be seen that the nucleon kinetic energies receive a large boost 
because of the high energy of the projectile. The angular distribution in the lab frame 
is obtained by integrating the photonuclear double differential cross section in the lab 
frame over all lab energies. In the projectile frame, the photonuclear angular distribution 
is isotropic; while in the lab frame it becomes non-isotropic and is peaked strongly in 
the forward direction. These features of the photonuclear angular distribution in the lab 
frame can be observed in figure 14. 

The nucleus - nucleus differential cross sections are obtained by taking the corre- 
sponding photonuclear cross section and integrating over the virtual photon spectrum, 
as discussed previously. It can be seen that all the nucleus - nucleus differential cross 
sections, in both the projectile and lab frames, follow the shapes of the corresponding 
photonuclear differential cross sections. This makes sense because the photonuclear dif- 
ferential cross sections are just integrated over the virtual photon spectrum. In the lab 
frame, the nucleus - nucleus differential cross sections, just like for the photonuclear dif- 
ferential cross sections, have nucleon kinetic energies that receive a large boost and angles 
that strongly peak in the forward direction. 


6 Comparison to experiment 

There is very little experimental data concerning differential cross sections for elec- 
tromagnetic dissociation. The best data available has been measured by Barrette et al. 
[39, 40], but much of their data involved spectral distributions of excitation energy. This 
will be analyzed in future work. However, a notable feature of their measurements is that 
all of their angular distributions are approximately isotropic in the projectile frame, which 
agrees with the assumption of the present work. Some kinetic energy distributions have 
been measured for outgoing neutrons and protons. See figure 13 of Barrette et al. [40]. 
In figures 21 and 22, we have compared our theory to these data. There are two points 
to note about these figures. Firstly, the experimental work quoted arbitrary units, so it 
was necessary to fit the absolute value (peak cross section) to the experiment. Secondly, 
the best fit is obtained by choosing the nuclear temperature constant in equation (13) as 
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d = 20. The comparison between theory and the limited experimental data is good. A 
better fit can probably be obtained with a more sophisticated approach to calculating the 
nuclear temperature and the spectral distribution (25). 


7 Conclusions 

This paper presents the first calculations of differential cross sections for electromag- 
netic dissociation in nucleus - nucleus collisions. These cross sections will be used in three 
- dimensional transport codes to improve estimates of radiation dose inside spacecraft. 
The results are given in a form in which they can be immediately used in fully three- 
dimensional space radiation transport codes that require differential cross sections in the 
lab frame. Cross sections are isotropic in the projectile frame and are in agreement with 
experiment. Spectral distributions in the projectile frame are compared to experimental 
results and are found to be in good agreement. However, the amount of data availble 
is very limited and the acquisition of more data would be very useful. Future work will 
involve a more sophisticated approach to calculating the spectral distribution and the 
nuclear temperature. 
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(A-l)o 


Figure 1: Photonuclear reaction with a two - body final state via the reaction 7 + A — > 
(A — l)o + N. This shows the decay of an excited parent nucleus, with nucleon emission, 
to the ground state of the daughter nucleus. The parent is excited to only one energy 
level E, which is determined by the incident photon of energy E 1 in the photonuclear 
reaction. The energy of the emitted nucleon is fixed. 



Figure 2: Photonuclear reaction with a three - body final state, 7 + A — » {A — 1) + 7 + N. 
This shows the decay of an excited parent nucleus, with nucleon emission, to a variety of 
energy levels in the excited daughter nucleus. The parent is excited to only one energy E, 
which is determined by the incident photon of energy E 1 in the photonuclear reaction. The 
energy level distribution p(E) in the daughter nucleus gives rise to a variety of energies 
E for the emitted nucleon, which gives rise to the Boltzmann distribution e~ E / k<s> . 
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Figure 3: Nucleus - nucleus reaction, A P + A T — > A P + (A T — 1) + N. This shows the 
decay of an excited parent target nucleus, with nucleon emission, to the ground state of 
the daughter nucleus. The parent can be excited to a variety of energies because the 
projectile nucleus A P brings in a variety of incident photon energies. The energy of the 
emitted nucleon can vary. 



Figure 4: Nucleus - nucleus reaction, A P + A T — > A P + ( A T — 1) + N . This shows the 
decay of an excited parent nucleus, with nucleon emission, to a variety of energy levels in 
the excited daughter nucleus. The parent can be excited to a variety of energies because 
the projectile nucleus brings in a variety of incident photon energies. The energy of the 
emitted nucleon can vary. 
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Figure 5: Relation between projectile and lab angles for an emitted photonucleon with 
kinetic energy T/v = mjv or a = 1.15259. The calculation used equation (67). 
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Figure 6: Relation between projectile and lab angles for an emitted photonucleon with 
kinetic energy T N = lOmjv or a = 1.00232. The calculation used equation (67). 
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Figure 7: Relation between projectile and lab angles for an emitted photonucleon with 
kinetic energy T N = 16m n or a = 0.999904. The calculation used equation (67). 
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Figure 8: Relation between projectile and lab angles for an emitted photonucleon with 
kinetic energy Tjv = lOOmjv or a = 0.998222. The calculation used equation (67). 
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Photonuclear (Projectile Frame) 



Figure 9: Photonuclear spectral distribution in the projectile frame evaluated at a photon 
energy of 20 MeV. The calculation was done using equation (14). 
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Figure 10: Photonuclear angular distribution in the projectile frame evaluated at a photon 
energy of 20 MeV. 
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Figure 11: Photonuclear double differential cross section in the projectile frame evaluated 
at a photon energy of 20 MeV. (Bottom figure is same as top figure, except rotated.) The 
calculation was done using equation (29). 
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Photonuclear (Lab Frame) 
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Figure 12: Photonuclear double differential cross section in the lab frame evaluated at 
a photon energy of 20 MeV. (Bottom figure is same as top figure, except rotated.) The 
calculation was done using equation (72). 
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Photonuclear (Lab Frame) 



Figure 13: Photonuclear spectral distribution in the lab frame evaluated a photon energy 
of 20 MeV. The calculation was done using equation (72) and integrating over the angle. 


Photonuclear (Lab Frame) 



Figure 14: Photonuclear angular distribution in the lab frame evaluated at a photon 
energy of 20 MeV. The calculation was done using equation (72) and integrating over 
energy. 


27 








AA (Projectile Frame) 



Figure 17: Nucleus - nucleus double differential cross section in the projectile frame. 
(Bottom figure is same as top figure, except rotated.) The calculation was done using 
equation (48). 
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AA (Lab Frame) 
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Figure 18: Nucleus - nucleus double differential cross section in the lab frame. (Bottom 
figure is same as top figure, except rotated.) The calculation was done using equation 
(48) and by transforming to the lab frame. 
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Figure 20: Nucle 
done using equat: 
energy. 




AA (Projectile Frame) 



Figure 21: Comparison between theory and experiment for proton kinetic energy distri- 
bution in the projectile frame. The reaction is 28 Si + Pb — » lp + 27 Al + Pb at 14.6 A GeV. 
Cross section units are arbitrary. Experimental data are from figure 13(b) of reference 
[40] and are represented by the large dots. Error bars are smaller than symbol sizes. 


AA (Projectile Frame) 



Figure 22: Comparison between theory and experiment for neutron kinetic energy distri- 
bution in the projectile frame. The reaction is 28 Si + Pb — » In + 27 Si + Pb at 14.6 A GeV. 
Experimental data are from figure 13(c) of reference [40] and are represented by the large 
dots. Cross section units are arbitrary. Error bars are smaller than symbol sizes. 
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